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ABSTRACT 



We present here a detailed spectral study of the X-ray emission of the persistent 
source and the low-fluence bursts of SGR J0501+4516 observed during a deep XMM- 
Newton observation near the peak of its 2008 outburst. For the persistent emission we 
employ a physically motivated continuum emission model and spectroscopically deter- 
mine important source properties; such as, the surface magnetic field strength and the 
magnetospheric scattering optical depth. We find that the magnetar surface temper- 
ature near the peak of its activity is 0.38 keV, corresponding to an emission area of 
131 km 2 at a distance of 2kpc. The surface magnetic field strength determined spectro- 
scopically, B = 2.2 x 10 G, is consistent with the dipole field strength inferred from 
the source spin and spin down rate. We fit the stacked spectra of 129 very faint bursts 
with a modified blackbody model and find a temperature of 1.16 keV, corresponding to 
an emission area of 93 km 2 . We also find an evidence for cooling during the burst decay 
phase. 
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1. Introduction 

Magnetars are manifestations of isolated neutron stars. Commonly known as Soft Gamma 
Repeaters (SGRs) and Anomalous X-ray Pulsars (AXPs), magnetars are generally characterized 
by relatively slow spin periods (in a narrow range of 2 — 12 s), rapid spin-down rates (5 x 10~ 13 — 
10~ 10 s s _1 ), and intense X-ray bursts with luminosities anywhere between 10 37 — 10 45 ergs -1 . Dipole 
magnetic field strengths of th ese sources, as inferred f rom t heir spin periods and spin-down rates, 
are m the order of ~ 10 14 G (IKouveliotou et alJll99fiL Il99flh . which are in perfect agreement with 



the predictions of the magnetar model that the re exist isolated neutron stars, which are powered 
by th eir own extremely strong magnetic fields (IDuncan fc Thompsonl [1992J; iThompson fc Duncan 
1993h . 
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Magnetars are persistent X-ray emitters in quiescence most of the time at luminosity levels 
around ~ 10 33 — 10 36 ergs -1 . As they enter active episodes, they usually exhibit significant varia- 
tions i n their persistent emis sion spectral and temporal behavior in conjunction with the outburst 
onset ( Rea fc EsDOsitdl201lh . Transient magnetars are generally undetectable unless they undergo 
outbursts, during which their X-ray fluxes are enhanced by factors of a few 10 — 100, followed b y 
a gradual decay over timescales of months to years (e.g. bogus et aliboid : iRea fc Espositolboilh . 
Typical magne tar bursts are short, lasting only for ~ 100 ms, with a peak luminosity lower than 



10 41 ergs _1 (jGogiis et al.ll200ll : Ivan der Horst et al.ll2012l ). Mor e energetic bursts occur much less 



frequently and have been obse r ved only from five SGR sources (jMazets et al.l 1 1973 ; iHurlev et al 



1999 



sec 



2005: 



Mazets et al 



1999 



Woods & Thompson! (12006^ and Mereghettil (120081 ). 



Watts et al 



2010; 



Gqgiis et al.ll201ll ). For reviews on magnetars, 



Both persistent X-ray emission and bursts are attributed to the decay of extreme l y strong mag- 
netic fields in the framework of the magnetar model ( Thompson &; DuncanlE"993 . 1995 : Thompson et al 
2002). However, it is not trivial to describe the observed data with detailed physical models, con- 
sidering that the radiation ought to emerge from the strong magnetic and gravitational fields. The 
persistent X-ray spectrum of magnetars is usually fitted with phen omenological models: either the 



combination of a blackbody and a power law (see secti on 14.3.1 of IWoods fe Thompson! I2006T). or 



two bl ackbody components (|Mereghettill2008l ) . Recently. I Giiver et al.l (|2007l . 120081 ) and lGogiig et al 
(|201ll ) studied the persistent emission spectra with an idealized physical model, the Surface Thermal 
Emission and Magnetospheric Scattering (STEMS). The STEMS mod el accounts for th e transport 
of radiation in a fully ionized and extre mely magnetized atmosphe re (jOzell l200ll . 120031 ) as well as 
for resonant cyclotron scattering (RCS) ( Lyutikov &: Gavriill 20od ) of photons emerging from the 
surface in the magnetosphere of the magnetar (see also Section [2]). A typical magnetar burst mainly 
emits radiation up to 100 keV. Its spectrum is well described phe nomenologically by the sum of tw o 
blackbody functions with temperatures o f ~3 and ~12 keV (e.g. llsrael et al.ll2008l : iLin et al.ll2012l ). 



To physically describe the burst spectra, iLvubarskvl (|2002l ) introduced analytically the shape of a 
thermalised spectrum modified by a strong magnetic field (~ 10 14 G). We use here the latter model 
in our burst spectral analysis. 

SGR J0501+4516 emitted a series of short X-ray/soft 7— ray bursts for almost two weeks start- 



and Suzaku 


— , 

Enoto et al. 


Lin et al. 


2011 


). A pointing 



the long term monitoring with RXTE and Swift /XRT , a spin-down rate of ~ 5.8 x 10~ 12 ss _1 was 
later obtained pro viding an estimate of the inferred dipole magnetic field strength of ~ 2 x 10 G 
(jGogiig et al.ll2010i ). Chandra observations of SGR J0501+4516 located the point source very pre- 
cisely at R.A. = 05 h 01 m 06!76, Dec = +45°16'33'.'92 (J2000), with a la uncertainty of 07 11 



(jGogiis et al.ll2010l ). Unlike many other galactic magnetars, SGR J050 1+4516 is locat ed in the anti- 
Galactic center direction, most likely at the Perseus arm at ~ 2kpc (|Xu et al.ll2006l ). As a result, 
the interstellar Hydrogen column density, Ah, towards SGR J0501+4516 is much less compared to 
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the values of magnetars located in the Galactic center direction, making SGR J0501+4516 an ideal 
source for studying its soft X-ray emission characteristics. 

In this paper we introduce for the first time the RCS process into the modified thermal burst 
emission, to study the burst spectra below lOkeV from SGR J0501+4516. We also analyze the 
persistent emission of SGR J0501+4516 with the STEMS model, and investigate whether there 
exists a link between the persistent emission and the dim burst spectra assuming the same scattering 
process. We describe the physical models used for the persistent emission and burst spectra in 
Section [2j In Section [31 we describe the analysis of the XMM-Newton data used here, and present 
the results of our spectral studies in Section HI Finally, we discuss the interpretation of these results 
in Section [5j 



Models for the Persistent and Burst Emission 



2.1. Persistent emission 



In the magnetar model, the quiescent X-ray emissio n could originate from the n eutron star 
surface heating by the decay of the strong magnetic field ([Thompson &; Duncanlll993l ). Persistent 
emission photons ought to go through the atmosphere and magnetosphere of the highly magnetized 
neutron star before they escape. The transport of radiation for different polarization modes was 
investigated in detail in a series of stu dies considering t he effects of absorption, emission , and 
scatt ering by the fully ionized plasma (jZane et al.1 boOll : |(Z)zell boOll ; lllo Lail boOll . bood ; lOzel 
2003). These showed that as photons propagate from the inner to the outer layers of a magnetized 
atmosphere, the spectrum of the emergent emission becomes harder than the original blackbody 
shape. 

As the phot ons reach the magnetos phere, resonant cyclotron scattering could effectively modify 
their spectrum ([Thompson et al.l 120021 ) . This process will again make the spect rum harder, by 
scattering the low energy photons to higher energies. iLvutikov &: Gavriill (|2006l ) calculated the 
simplified one dimensional resonant cyclotron scattering of a Planck spectrum by a non-relativistic 
warm plasma in an inhomoge neous magnet i c field. Ap plying the same scattering process to the 
emission from the atmosphere, iGiiver et al.l (|2007l . l2008l ) developed the Surface Thermal Emission 
and Magnetospheric Scattering model (STEMS). They then used this model to successfully describe 
the persistent emission spectra of XTE J1810— 197 and 4U 0142+61, and also estimate the surface 
magnetic fields of these two magnetars. The STEMS model depends on four parameters: the surface 
magnetic field (B), the temperature of the neutron star (kT), the temperature of the plasma in 
the magnetosphere (/3), and the optical depth for resonant cyclotron scattering (r). Because of the 
strong gravitational field of the neutron star, the STEMS model also includes a fixed gravitational 
redshift (z) parameter, which depends only on the ratio of mass over radius of a neutron star. We 
fixed z = 0.306 throughout this paper, corresponding to a typical neutron star with mass of 1.4 Mq 
and radius of 10 km. 
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2.2. Burst emission 



One of the mechanisms invoked for the origin of the short duration magnet ar bursts is the sud- 
den r elease of energy by cracking of the magnetically strained neutron star crust ([Thompson Duncan 
19951 ). This sudden energy release creates a hot pair plas ma, which then becomes tra pped in the 
magnetosphere forming a bubble with large optical depth ([Thompson Duncanlll995l ). The pho- 
tons need to go through multiple scatt erings to escape the bubble, thereby thermalizing the burst 
radiation ([Thompson &: Duncan! Il995l ). Unlike in the case of the persistent surface emission for 
which the absorption, e mission and scatt ering are important, scattering is the most dominant pro- 
cess in the burst bubble (|Lyubarskyll2002l ) . Taking into account the fact th at photons with different 
energies would have different scattering cross-sections, iLvubarskvl (|2002l ) calculated the radiation 
transfer for the magnetar bursts. The resulting spectral shape has a much flatter low energy distri- 
bution compared with the equilibrium blackbody spectrum (hereafter we call the resulting shape as 
the modified blackbody). The modified blackbody photons are also subject to multiple scatterings 
as they go through the magnetosphere. We a dopt ed the same proces s of the resonant cyclotron 
scattering as used in iLyutikov &; Gavriill (|2006l ) and iGiiver et al.l (120071 ) to up-scatter the modified 
blackbody burst spectrum. The scattered modified blackbody and the Planckian function differ sig- 
nificantly at energies below the bolometric temperature (T&). The softer emission is highly affected 
by the interstellar absorption, which would make the two models indistinguishable if the source 
were towards a high interstellar iVn region. The location of SGR J0501+4516 makes it, therefore, 
an excellent source to determine the nature of its burst emission mechanism. We generated a 
numerical grid for the modified blackbod y with the res onant cyclotron scattering (MBB+RCS) 
to fit the burst spectral data in XSPEC (|Arnaudl ll996), The MBB+RCS model has three free 
parameters: the bolometric temperature of the bubble (kTf,) which was allowed to range between 
0.1 — 20keV, and f3 and r as defined in Section 2.1. We also used the same gravitational redshift 
parameter (z = 0.306) while fitting the burst spectra. 



3. Observations and Data Reduction 



SGR J0501+4516 was observed with the XMM-Newton Observatory (jjansen et al.l l200ll ) in 
seven occasions between August and September 2008. Here, we selected only the first observation 
performed on 2008 August 23 (Observation ID: 0560191501), the most burst active day of the 
source, to study the spectral properties of the bursts and of the underlying persistent emission 
at the same time. The observation lasting for 48.9 ks was perform ed in the small win dow mode 
of the European Photon Imaging Camera (EPIC) pn instrument (jStriider et al.ll200ll ). allowing 
a temporal resolution of 6 ms. We processed the data using SAS version 11.0.0 with the latest 
calibration files generated on 2012 May 18. 

We identified and removed events in the piled-up time intervals (64 s in total), which are 
the ones with count rates over 50 counts s -1 , from the source event list. We then constructed the 



- 5 - 



source lightcurve with 100 ms binsize as shown in Figure 1. The average count rate of the persistent 
emission is ~ 5.5 counts s" 1 , corresponding to no more than 5.5 counts per bin in the light curve of 
Figure 1 (or a count rate of 55 counts s" 1 ). We assumed that a bin contained a burst if it exceeded 
the persistent X-ray emission level by at least 2a, corresponding to ~ 10 counts per bin (or a count 
rate of 100 counts s _1 ). We, therefore, accumulated the persistent emission spectrum using 100 ms 
time bins with count rate less than 50 counts s -1 , and the burst spectra from those bins with over 
100 counts s -1 . The total exposure time of the persistent emission spectrum is 32.7 ks. There are 
129 identified bursts having 100 ms peak count rates between 100 and 500 counts s _1 ; the total time 
encompassed by these bursts was 8.7 s. We combined all 129 burst data into one stacked spectrum 
and regrouped the resulting spectra to have a minimum of 15 counts in each spectral bin, making 
sure that we kept the bin size within the spectral resolution of the instrument. We then fit the 
spectra using XSPEC version 12.7.0. 
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Fig. 1. — The lightcurve of SGR J0501+4516 in 100 ms temporal resolution. The dashed line marks 
the count rate (50 counts s _1 ) selected for the persistent emission level. The dotted line denotes 
the count rate threshold for the burst selection. 
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4. Spectral Analysis Results 
4.1. The Persistent Emission Spectrum 



We first fit the spectrum of the persistent emission with a single blackbody and a power law 
(BB+PL) . The best fit parameter s liste d in Table [H are consistent with the results reported in 



Rea et all (|2009l ) and lGogiis et al.1 (|2010h . We then fit the spectrum with STEMS (see also Figure 



|2J). The STEMS model can fit the persistent spectrum equally well as the BB+PL model, and 
provides well constrained model parameters (see Table [1]). We find the unabsorbed source flux to 
be (5.88 ±0.02) x 10~ n ergs" 1 cm" 2 in the 0.5-6.5keV range, and a temperature kT = 0.38 ±0.02 
keV, corresponding to a hot spot surface area of 131 ± 27 km 2 (assuming a source distance of 2 kpc). 
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Fig. 2. — The count spectrum of the persistent emission of SGR J0501±4516 overplotted with the 
best fit STEMS model curve (upper panel), and the fit residuals (lower panel). 



21 



Using the STEMS fit, we obtained a hydrogen column density of = (6.7 ±0.2) x 10" cm 
that is much lower than the one obtained with the BB±PL fit (Nn = 9.1 x 10 21 cm~ 2 ). The 
weighted average va lue of iNfa towards th e direc tional of SGRJ0501+4516 is 6.2 x 10 21 cm~ 2 and 5.2 x 
10 21 cm -2 using the lDickey fc Lockmanl (1990) and the lKalberla et aJj ( 20051 ) surveys, respectively. 
The Njj from the STEMS fit agrees well with the Galactic value, however, the BB±PL model fit 



x As obtained with the HEASARC nH tool: http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 



-7- 



overestimates the absorption. 



4.2. The Time-integrated Burst spectrum 



To remove the contribution of the persistent emission from the burst emission, we used the 
persistent spectrum as the background for the burst spectrum. We fit the burst spectrum with sev- 
eral continuum models: blackbody, modified blackbody, blackbody with RCS, and MBB+RCS. In 
all our fits, we fixed N# at 6.7 x 10 21 cm -2 , the value obtained from the persistent spectrum fit with 
the STEMS model. Previous studies revealed that the magnetospheric paramet ers do not vary sig- 



nificantly over timescales of months (|G fiver et al 



2007 



20081 ; iGogiis , et al.ll201ll ). We also assumed 



that the RCS by the non-relativistic warm plasma in the magnetosphere of SGR J0501+4516 is not 
varying over the course of the XMM-Newton observation and adopted the magnetospheric electron 
velocity (/?) and the optical depth (r) from the persistent spectrum fit. We list all the fit results 
in Table [TJ Figure [3] exhibits the burst spectrum with the MBB+RCS fit (top panel) and the 
residuals for all models (next four panels). We find that a MBB+RCS gives the best fit among all 
four models. Note the apparent deviation above 7keV in all spe ctral fits; t his m ay be due to an 
additional component in the spectra of SGR bursts as shown by lLin et al.1 (J2012J). We calculated 



the average unabsorbed flux of the bursts in the 0.5 — 10 keV to be (1.80 + 0.05) x 10 9 ergs" 
The total energy released by all 129 very dim bursts is ~7.5xl0 36 erg. 
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4.3. The Time-resolved Burst spectrum 

We studied the time evolution of the spectral properties of the brightest bursts in our set of 
129 events, using their stacked time-resolved spectra. We selected all the bursts with peak count 
rate more than 300 counts s _1 on a 25 ms binsize. Note the fact that we reduced the time resolution 
to accommodate SGR bursts with finer temporal structures. This selection resulted in a set of 47 
bursts. Then, we separated each selected burst into three parts: rise, peak, and decay. The peak 
bin is defined as the 25 ms time bin with the highest count rate during the burst; the rise part 
is the time interval from the onset of the burst to the start of the peak bin; and the decay part 
starts from the end of the peak bin ending when the emission returns back to the background level. 
Figure U] shows an example of the three parts of a relatively bright burst. We generated stacked 
spectra for the rise, peak and decay parts and fit these spectra with the MBB+RCS model. 

Figure [5] exhibits the spectra of the rise, peak and decay parts of the bursts, and their best fit 
model curves. We list all spectral fit results in Tabled) We find that the burst bubble temperature 
is the highest at the peak among the three parts. However, the temperatures of the rise and 
peak parts differ only by ~ 2a. Therefore, we fit these two parts simultaneously, forcing them 
to yield a single (best fit) temperature of 1.25 ± 0.08 keV. We notice that the temperature of the 
modified blackbody is lower (> 3a) during the decay phase of the bursts. The average unabsorbed 
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Fig. 3. — The time-integrated stacked spectrum of 129 dim SGR J0501+4516 bursts overploted 
with the best fit MBB+RCS model and the residuals for MBB+RCS, BB+RCS, MBB, and BB 
models (from the top to bottom). 



burst flux (0.5 - lOkeV) is (9.8 ± 0.7) x 10" 10 ergs -1 cm" 2 , (5.7 ± 0.4) x 10" 9 ergs" 1 cm" 2 , and 
(9.1 + 0.5) x 10" 10 ergs" 1 cm -2 for the rise, peak and decay parts, respectively. We present a coarse 
time evolution of the burst temperature and flux in Figure 



5. Discussion 

We investigated the persistent X-ray and very dim burst emission properties of SGR J0501+4516 during 
the day after the onset of its 2008 outburst, using the deepest XMM-Newton observations of 
the source that included a large number of bursts. We analyzed the persistent X-ray emission 
with the STEMS magnetar emission model and found a magnetar average surface temperature 
of kT = 0.38 keV, that is much lower than the temperature of ~ 1 keV, obtained with the phe- 
nomenological BB+PL model fit. The persistent emission of almost all other sources modeled 
with BB+PL yields a blackbody temperature of ~0.5keV. Note that this value is much lower 
than that of SGR J0501+4516 fit with BB+PL, but closer to or slightly higher than that esti- 
mated with the STEMS model fit. We also found that for SGR J0501+4516, the BB+PL model 
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Fig. 4. — An example of a bright burst lightcurve from SGR J0501+4516 with 25 ms time resolution. 
The black bar is the peak bin. The rise and decay parts are filled with dots and horizontal 
histograms, respectively. The reference time of the x-axis is 2008 August 23 12:19:16.861 UTC. 



10^ cm 2 than the Galactic value (5— 6x10 cm 



gives a larger Nn ~ ]n22 "•" 1 
SGR J0501+4516, while the STEMS N R 



-'- ~ >N i j n the direction towards 
6.7 x 10 21 cnr 2 agrees well with the latter value. 



Enoto et al.l (|2010i ) investigated the Suzaku observations of SGR J0501+4516 on 2008 August 
26 — 27 in the energy r ange of 1 — 70 keV. They u sed a thermal emission with a resonant cyclotron 
scattering model (see iLvutikov fc Gavriill 120061 ) to describe the persistent soft X-ray spectrum 
plus a simple power law for the hard X-ray component. This model did not provide a good fit 
(reduced \ 2 = 1-73). They reported a value for JVh = 6.2 x 10 21 cm -2 , a surface temperature 
of kT = 0.39 keV, and magnetospheric scattering parameters, r = 5.1, and /3 = 0.3, which are 
quite similar to our results. However, the uniqueness of the STEMS model is the fact that it 
takes into account the presence of the strong surface magnetic field and its possible effect on the 
emergent spectrum via atmospheric processes. Using STEMS we obtain, for the first time, the 
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Fig. 5. — The count spectra and their best fit with the MBB+RCS model curves for the rise (dots), 
peak (circles), and decay (triangles) parts of the stacked bursts. 



surface magnetic field strength of SGR J0501+4516 via X-ray spectroscopy, as B = 2.2 x10 G, a 
value which agrees well with the inferred equatorial dipole field strength of 2xl0 14 G (|Rea et al 
20091 : IGogiis et allhoioh . 



The STEMS model has already been used to study the persistent emission sp ectra of four 



AXP s , XTEJ1810-197, 4U 0142+61, IE 1048.1 - 5937 and 1RXS J170849.0- 400910 (IGiiver et al 



and SGR J1550 - 5418 (IGogiis et al 



20071 . 120081: lOzel et al.ll2008h. and four SGRs, SGR 1900 + 1 4, SGR 0418 + 5729, SGR 0526 + 66 



2011 



Giiver et al. 



2011 



2012 



Ng et al.ll201ll ). In most cases, 



including SGR J0501+4516, the surface magnetic field obtained from X-ray spectroscopy is con- 
sistent with the dipole magnetic field inferred from the spin period and the spin down rate, as 
shown in Figure [71 Because magnetars emit energetic bursts, their surface magnetic field strengths 
are expected to be in the 10 14 — 10 15 G domain . Some sources, however, may have lower inferred 
dipole magnetic field strengths (|Rea et al.ll2010l ). although higher order magnetic structures (mul- 
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Fig. 6. — Time evolution of the temperature of the modified blackbody (squares) and the unab- 
sorbed flux (filled circles) in the 0.5 — lOkeV range for the rise, peak and decay parts of the 47 
brightest bursts of SGR J0501+4516. The error bars in the time axis represent the total exposure 
time of each accumulated part. 



tipolar fields) may be present ( Giiver et al.ll201~ll ). The similarity between the surface magnetic 
field strengths of SGR J0501+4516 as obtained via X-ray spectroscopy and as inferred from the 
spin properties suggests that its surface magnetic field topology is close to a dipole. 

We also studied here, for the first time, the dimmest bursts and their physical link to the 
persistent X-ray emission as both kinds of emission are redistributed by resonant cyclotron scat- 
tering process in the magnetosphere. The burst spectrum can be described with a warm bubble 
trapped by the magnetosphere, emitting therm al radiation modi fied by a strong magnetic field. 
The bolometric temperature Tf, is introduced by iLvubarskvl (|2002l ) with the total radiation flux as 



F = 2 a Th- The emission area of the burst can be approximately calculated as 



A 



A-kD 2 f 

2 al b 



(1) 



where A is the emission area, D is the distance to SGR J0501+4516, / is the unabsorbed flux, a 
is the Stefan-Boltzmann constant, and T& is the bolometric temperature. Assuming a distance to 
the source of ~ 2kpc, we estimate the average emission area of the burst as 93 ± 10 km 2 , which is 
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Fig. 7. — The comparison of the inferred dipole magnetic fields with the surface magnetic field 
obtained with X-ray spectroscopy for nine magnetars. The dashed line marks the line on which the 
two axes are e qual. SGR J0 5 01+4 516 is shown as a five-point star. All other values in this figure 
are taken from iGiiver et al.l (|201ll ). SGRJ0418 + 5729 has an upper limit for the inferred dipole 
magnetic field and a measurement with STEMS spectroscopy. 



7.4% of the magnetar surface (Rns = 10 km). 



Enoto et al.1 (|2010l ) and iNakagawa et al.l (|201ll ) studied the spectra of the persistent emission 
and the dim short bursts (average flux ~ 10 -9 ergs -1 cm -2 ) from SGR J0501+4516 with Suzaku. 
They fit both spectra with two blackbqdy f unctions and a power law model. Comparing the 
parameters of these fits. INakagawa et al.l (|201ll ) concluded that the persistent emission and the dim 
bursts have the same emission mechanisms. They further suggested that the persistent emission 
of magnetars may be composed of numerous microbursts. We find that both kinds of emission 
are forms of modified thermal emission and redistributed by the resonant cyclotron scattering in 
the magnetosphere. However, the spectral shapes of the thermally dominated portion (<4 keV) of 
the STEMS and MBB+RCS models are remarkably different and they yield significantly different 
temperatures. We, therefore, argue against the idea that the persistent and burst emissions have 
a common origin. The main difference between the two is that the scattering is the dominant 
process in the burst bubble, while for the persistent emission, one needs to consider not only the 
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scattering but also the emission and absorption processes as photons travel through the magnetar 
atmosphere. We investigated the dependence of burst peak times on the rotational phases of the 
source and we found that bursts occur quite uniformly over the entire spin phase, showing no 
significant concentration during any phases. This result supports our argument that short bursts 
and persistent emission do not have a common origin. 

Finally we investigated the stacked time-resolved spectra of the brighter bursts from SGR J0501+4516. 
We found that their temperature traces their flux behavior and has a lower value in the decaying 
part compared t o the rise and pe ak parts. We note here that these results do not contradict the 



earlier results bv lLin et al.l (|201ll ). who find the burst hardness to evolve from hard-soft-hard with 
flux over a very broad burst flux range, since the bursts studied here correspond to a very narrow 
flux range in a softer energy band. From the Equation 1 above, we calculate the emission area for 
the rise, peak and decay parts as 58 ± 19 km 2 , 160 ± 41km 2 , and 118 db 29 km 2 , respectively. The 
emission area during the peak of the burst is similar to the area of the persistent emission, within 
lcr. Although more coarse, the behavior of cooling during th e decaying tail is s imilar to what has 



been seen in the tai ls of energetic bursts from SGR 1900 + 14 (jLenters et al.ll2003l ) and SGR 1806-20 



(jGogiig et al.ll201ll ). implying that cooling takes place at all energy scales of SGR events. 



We thank Yuri Lyubarsky for stimulating discussions. L.L. is funded through the Post-Doctoral 
Research Fellowship of the Turkish Academy of Sciences (TUBA). 
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Table 1. Spectral fit results of the persistent and burst emission from SGR J0501+4516. 



Spectrum Model N H B kT & index r xt/dof 

(10 22 cm- 2 ) (10 14 G) (keV) 



Persistent BB+PL 0.91 ± 0.01 - 0.70 ± 0.01 2.79 ± 0.04 - - 0.7657/117 

Persistent STEMS 0.67 ±0.02 2.21 ± 0.07 0.38 ± 0.02 - 0.37 ± 0.01 5.0 ± 0.2 0.7615/116 



Burst all 


BB 


0.67 


1.08 + 0.02 






2.1054/84 


Burst all 


MBB 


0.67 


1.35 ± 0.04 






1.4944/84 


Burst all 


BB+RCS 


0.67 


0.93 ± 0.03 


0.37 


5.0 


1.4415/84 


Burst all 


MBB+RCS 


0.67 


1.16 + 0.04 


0.37 


5.0 


1.174/84 


Burst rise 


MBB+RCS 


0.67 


1.12 + 0.11 


0.37 


5.0 


1.019/20 


Burst peak 


MBB+RCS 


0.67 


i oc+0.12 
±.oo_ .i 


0.37 


5.0 


0.9632/25 


Burst decay 


MBB+RCS 


0.67 


92+ 06 
— 0.07 


0.37 


5.0 


1.169/36 


Burst rise+peak MBB+RCS 


0.67 


i or+0.09 

L - zo -0.07 


0.37 


5.0 


1.015/46 



a For all burst fits this parameter is the feTj, of the modified blackbody model 



